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Stability about cosmological background solutions to the bi-metric Hassan-Rosen theory is studied. 
The results of this analysis are presented, and it is shown that a large class of cosmological back- 
grounds is classically unstable. This sets serious doubts on the physical viability of the Hassan- Rosen 
theory — and in turn also of the de Rham-Gadabaze-Tolley model, to which the mentioned theory 
is parent. A way to overcome this instability by means of curvature-type deformations is discussed. 
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Introduction & Overview — The general theory of rel- 
ativity proposed by Einstein in 1916 [1] provides the fun- 
damental building block of our current understanding of 
gravitation. This framework — describing the dynamics 
of a massless spin- two field in four dimensions — has been 
tested from scales of about a fraction of a millimeter up 
to scales of a few astronomical units, and agrees remark- 
ably well with all experimental data. 

Despite of its successes, and the necessity of a theory 
of quantum gravity in the ultra-violet, it remains rather 
unclear whether general relativity is a valid description 
also on cosmological scales. Therefore it is tempting to 
study its consistent infra-red deformations. Several of 
those possibilities have been considered such as extra- 
dimensional models [MS], multi-gravitation [MS], or de- 
formed (e.g. massive) gravity [5] I9HI3]. 

Since the fundamental work of Fierz and Pauli [TT] 
in 1939, who constructed a consistent theory of mas- 
sive gravity on Minkowski background to linear order, 
the quest has long been unsuccessful to consistently gen- 
eralize such a framework to curved space-times. In 
Ref. [13] this task has been established on a Friedmann- 
Lemaitre-Robertson- Walker background, which — by in- 
clusion of the Ricci scalar — was shown to be full re- 
spected throughout the entire realistic cosmological evo- 
lution. An important feature of this theory is that the 
Fierz-Pauli mass parameter can be consistently set to 
zero, therefore providing a modification of general rela- 
tivity solely on curved space-times. This might be very 
important in the light of the Boulware-Deser ghost [14] 
and the vdvz discontinuity [131 HI] ■ 

Many of the models which have been proposed so far 
to modify gravity have the unphysical need to fix a refer- 
ence metric, or, if this metric is dynamical, lack the exis- 
tence of a respected cosmological background. A recent 
attempt to modify gravity with a bi-metric theory which 
allows for a cosmological background has been presented 
in [8]. This work was only concerned with establish- 
ing a realistic cosmological background. A complete and 
consistent study of fluctuations about this background is 
very important for stability issues (c.f. [TT]). 

In this work we present results of precisely such a sta- 
bility analysis, and show that a large class of the cos- 



mological branch of the Hassan-Rosen theory [7] is not 
physically viable. We then show (for one particular case) 
a way to ensure full stability, at least on the linear level. 
This route, which has been already paved in Ref. [IB"] , 
might be key towards a viable non-linear completion. 

Framework — The bi-metric action under consideration 
is (c.f. Ref. 0) 
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Here, X := y/ g~ 1 f, M. is a four-dimensional pseudo- 
Riemannian manifold, the metrics / and g have signature 

( — , +, +, +), and the units are such that % = c = 1. For 
the sake of convenience, the matter fields — which are 
minimally coupled to g in the matter Lagrangian C m — 
are denoted by $. R[ • ] is the Ricci scalar of the respec- 
tive metric, the j3 n are fixed, real parameters, and e n (X) 
are elementary symmetric polynomials of the eigenvalues 
of the matrix X, e.g. 



eo(X) = 1 , 
e 2 (X) = i(fl 
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where the double-lined square brackets denote the matrix 
trace, i.e. [X] = X M M . The quantity e^ is not displayed as 
it will not be included in the present analysis, which refers 
to the cosmological studies of Ref. [5] wherein /3 3 is set 
to zero. Actually, in the case of massive gravity (where 
/ is non-dynamical) one can show that its inclusion is 
phenomenologically non-acceptable [TH]. Hence, we will 
set /?3 = 0. It is easy to check that the model (except 
the matter sector) only depends upon three dimensionless 
parameters [2D]: 
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Stability — To check for stability or instability, respec- 
tively, one has to expand the fields / and g about certain 
backgrounds which are consistent with the action (JlJ. 
Then one studies how the perturbations evolve. 

To this end we expand / and g about the backgrounds 
/(o) and 3(0) j respectively, 



/ = /(o) + Sf 
9 = 5(o) + Sg . 



and define the matrix 8 via 
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(4b) 
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which occurs in the whole interaction term in ([!]), and 
allows to express /m) through g( ) y i & /(o) — 9(0) ■ 
In the cosmologically-relevant case the background 
of the fluctuation dg (to which our matter sector is 
coupled to) has a Friedmann-Lemaitre- Robertson- Walker 
structure, i.e. gm-\ = diag(— 1, a 2 , a 2 , a 2 ), where a is the 
scale factor, being normalized such that it equals one to- 
day. Then one has 



= diag(|a|,|/3|,|/3|,|/3|), 



(6) 



wherein a = a(a) and ft = ft{a) depend on a. In general, 
both backgrounds are not the same, i.e. a (a) ^ 1 and 
ft (a) ^ 1, and one has (c.f. Ref. [5]) 



a(a) 



_ d(aft(a)) 
da 
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The function /3(a) is given by a root of a quartic poly- 
nominal, determined by the Friedmann equations. How- 
ever, for the choice of fts = 0, this equation is only cubic 
in ft. The corresponding three solutions show quite dif- 
ferent behavior. In fact, in the limit of small scale factor 
two of them diverge, while one approaches zero. The lat- 
ter solution implies that those terms, that are contracted 
with the inverse metric , have prefactors that diverge 
(in the limit of vanishing scale factor) the stronger the 
higher their order is. 

Actually, as wc will see below, fluctuations become of 
order one already at some moderately small value of a. 
Moreover, the two other solutions imply a zero-crossing 
of a(a), which makes 9{a) non-analytic — a particular- 
ity which also concerns the associated solutions for the 
fluctuations. 

Let us now come to the general cosmological 
discussed in Ref. [8]. This case is, in particular, deter- 
mined by the solution to ft corresponding to ft (a = 1) = 1 
|21j . which implies that the density parameter £1 := p/p c , 
equals one today, i.e. £l(a — 1) = 1, being demanded by 
CMB observations [T^. Here, p is the energy density of 



the Universe and p c is its critical value for which the 
spatial geometry is flat. 

By performing distance-related tests using cosmolog- 
ical data, it has been shown by the authors of Ref. [8] 
that it is possible to choose the above parameters such 
that a realistic cosmological background can be obtained. 
Unfortunately, this has been done only for a very limited 
range of red shifts, and is purely on the background level. 

Linearization of the action ([I]) is straightforward but 
yields a rather lengthy expression for the action of the 
fluctuations Sf and 8g. Its detailed form is not very en- 
lightening. Here, instead, we will schematically present 
the field equations and their full solution for certain com- 
ponents of the fluctuation — in comparison to the unde- 
formed theory, i.e. purely general relativity coupled to 
appropriate matter, plus one isolated tensor field. 

For the sake of studying stability of the backgrounds it 
will be sufficient to determine the behavior of the spatial 
components (ij) of the zero-mode fluctuations. For i =/= j, 
one finds schematically the set of coupled field equations 

6f"j + at Sf-j + b x Sfij = ci Sgij , (8a) 

Sg'-j + a 2 Sglj + b 2 Sg i2 = c 2 Sfij , (8b) 

wherein a prime denotes a derivative with regard to the 



scale factor a, and the quantities a^, bj, 
besides on a 
Defining 



depend — 

also on other parameters of the model. 
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one can express the system ( 8a b) as 



IJ Ulj 
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Stability of the above system depends upon the behavior 
of the real parts of the eigenvalues of the matrix A^- . 

Analyzing precisely those real parts of the mentioned 
eigenvalues shows, in particular, two things: First, and 
as a check, we confirm that undeformed theory is stable. 
Second, if the deformation is switched on, there are real 
parts of the eigenvalues of the matrix Ay which show 
a quite peculiar behavior: One diverges at some finite 
value a = a*, where the divergence is such that it grows 
towards +00 for (a — a*) — > + , and tends towards —00 in 
the case (a — a*) — > CU. Furthermore — again, contrary 
to the undeformed theory, which stays finite — there is 
an eigenvalue whose real part grows strongly for small a. 
More precisely, for a — > it diverges much faster than 
1/a, yielding an exponential divergence of the associated 
mode. 
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Figure 1: Absolute values of the relative metric fluctua- 
tions := 5f 23 /f^ > (red curve) and y% := 8g 23 /g^ (blue 
curve) as functions of the scale factor a. Dashed, colored 
lines correspond to the undeformed theory, solid lines to the 
deformed one. The dot-dashed, vertical line is at a = a* ~ 0.3 
(for the parameter set |T2|)). 



We should stress that all parameters within the physi- 
cally-relevant intervals (as given in Ref. [5]) 



1-5 < M* < 3.0 , 2 < M < 3.5 
-0.5 < /3 2 < -0.1 , 



(11) 



yield the same quantitative behavior, and that in all those 
cases the scale factor at which the theory is non-analytic, 
a*, is far larger than its value at recombination. On the 
other hand, it is smaller than the value up to which super- 
novae data have been analyzed in [8]. Exactly the same 
holds true for the value at which fluctuations become of 
order one, a n i (see below). Choosing 



M* = 2.5 , M = 3.0 , 



-0.3 



(12) 



we find a n i « 0.1 and a* « 0.3. 

Fig. [I] shows the full solution to the system (8a b). It is 
shown (exemplary for the parameter set ( 12 ), and certain 
initial conditions) how the (absolute values of the) rela- 
tive fluctuations behave as a function of the scale factor 
a. One can read off the aforementioned instability from 
the solid lines, describing one particular realization of 
the deformed theory. In contrast, the undeformed theory 
is well-behaved. One also observes the same unphysical 
backward instability as in Refs. |12j . implying that the 
set-up is merely self-protected. We checked that the in- 
stability occurs for all cosmologically-allowed parameters 



out of the intervals ( 11 ) (for the present case of no mat- 
ter coupling to /, as discussed in Ref. [H]). On top of 
that, it is independent of the precise details of the initial 
conditions. 



Figure 2: Stability parameter plot in the /3-j plane, exem- 
plary for a — 0, f3 2 — —0.4 and M = 1. Green dots represent 
fully stable regions, yellow ones indicate classical instability, 
and red points stand for unitarity violation (c.f. main text). 



Let us finally come to one particularly interesting case, 
which is constituted by the limit M* = Mf/M g — > oo. 
Therefor the /-field is frozen into its background value, 
which may be taken to be Minkowskian due to the lack 
of respective matter couplings (c.f. Eq. ([l}). Performing 
analogous studies as above reveals the same mentioned 
backward instability — the figure corresponding to Fig. [l] 
looks qualitatively the same in this respect. 

Since, now, there is only one dynamical metric (albeit 
with a particular deformation term) we can easily use 
a modified version of the stability analysis performed in 
Ref. [13] . This tantamounts to study — after introduc- 
tion of Stiickelberg fields — the roots of the determinant 
of the full kinetic operator, from which bounds for sta- 
bility and unitarity can be directly read off (c.f. Ref. [13] 
for details). Moreover, Ref. [13] also provides a guide- 
line towards a resolution to the mentioned instability is- 
sues. This way constitutes of adding to the Lagrangian 
curvature-type deformations of the form 



0]c 



H a 9q 



(13) 



which are unique up to second adiabatic order [13] . Here, 
a subscript indicates a go-t>ackground quantity, a, (3, 7 
are real dimensionless parameters, and square brackets 



4 



around indices stand for anti-symmetrization. Including 
terms of higher adiabatic order requires introducing fur- 
ther parameters with appropriate inverse mass dimension 
to compensate for the additional derivatives acting on g . 

Fig. [5] shows a parameter plot (/3-7 plane) for a = 0, 
/?2 = —0.4 and M = 1. One observes, in particular, two 
things: First, the case ofa = /3 = 7 = (which cor- 
responds to the original model) is classically unstable, 
albeit it does not violate unitarity (as expected), and 
provides an independent confirmation of the aforemen- 
tioned instability. Second, there exists parameter values 
(of order one) such that the theory is truly stable! This 
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inevitably involves curvature extensions (c.f. Eq. (13)) 



Although in our last analysis we only checked the case 
of frozen /, we expect similar results for the full model, 
in which both metrics are dynamical. We will devote a 
future publication to such an analysis in extended form. 

Summary & Outlook — Let us summarize: Starting 
form the general bi- metric model 0, using the gen- 
eral phenomenologically- viable parameters intervals (11) 
(given in Ref. [8 ) that allow for cosmological back- 
grounds, and expanding the action to second order in 
the fluctuations about these backgrounds, one finds that 
the theory under consideration is classically unstable. 

We confirmed (for a special case) our results with an 
independent analysis method (introduced in Ref. |13j). 
and were able to show — by appropriate supplementation 
with the curvature-type deformation terms (13) — that 



the theory can be made stable on the linear level. 

We expect that this is also possible for the case in 
which both metrics are dynamical. This should in par- 
ticular hold for the case of common cosmological back- 
grounds (i.e. when / is also coupled to matter). In fact, 
for the so-called minimal model — in which the combina- 
tion Sf — Sg acquires a Fierz-Pauli-like mass term — this 
seems clear (c.f. Ref. [13]). Because the Hassan-Rosen 
theory is parent to the de Rham-Gabadadze-Tolley model 
of massive gravity (which is obtained from the theory ([!]) 
in a certain limit), the physical relevance of massive grav- 
ity of this type seems also questioned. 

In order to cure the aforementioned instabilities, one 
has to think about extensions to the bi-metric model ([l} . 
As mentioned above, a prime starting point could be the 
model proposed in |13j which includes curvature terms to 
ensure classical stability over the entire cosmic expansion 
history, thus resolving the issue of self-protection (c.f. [12j 
113]). Although the mentioned model is only a one- metric 
theory, it might point in the right direction, and its bi- 
metric non-linear completion (which necessarily involves 
higher-derivatives of at least one of the spin-two fields) 
seems tempting to explore, albeit it might very well be 
that nature prefers massless gravity. 
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